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The cosmological model of dark energy interacting with cold dark matter without coupling to 
the baryonic matter, is studied in the background of both classical Einstein and loop quantum 
cosmology. We consider two types of interacting models. In the former model, the interaction is a 
linear combination of the densities of two dark sectors, while in the latter model, the interaction 
with a constant transfer rate depends only on the density of cold dark matter. It is shown that the 
dynamical results in loop quantum cosmology are different from those in classical Einstein cosmology 
for both two kinds of interacting models. Moreover, the form of the interaction affects significantly 
, the dynamical results in both kinds of cosmology. 
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I. INTRODUCTION 



Recently, the discovery of the acceleration of cosmological expansion at present epoch has been the most principal 
achievement of observational cosmology. Numerous cosmological observations, such as Type la Supernovae (SNIa) 
Cosmic Microwave Background Radiation (CMBR) |2| and Large Scale Structure [3j, strongly suggest that the universe 
is spatially flat with about 4% ordinary baryonic matter, 20% dark matter and 76% dark energy. The accelerated 
expansion of the present universe is attributed to the dominant component of the universe, dark energy, which has 
a large negative pressure but not cluster. In fact, it has not been detected directly and there is no justification for 
assuming that dark energy resembles known forms of matter or energy. A large body of recent work has focussed on 
understanding the nature of dark energy. However, the physical origin of dark energy as well as its nature remain 
enigmatic at present. 

The simplest model of dark energy is the cosmological constant A [J] , whose energy density remains constant with 
'— time pa = A/8irG (natural units c = H = 1 is used throughout the paper) and whose equation of state (defined 
as the ratio of pressure to energy density) remains w = — 1 as the universe evolves. Unfortunately, the model is 
burdened with the well-known cosmological constant problems, namely the fine-tuning problem: why is the energy 
of the vacuum so much smaller than its estimation? and the cosmic coincidence problem: why is the dark energy 
\ density approximately equal to the matter density today? These problems have led many researchers to try different 
C*") ■ approaches to the dark energy issue. A possible method is to assume the equation of state (EoS) w is a dynamical 
variable, and thus the dynamical scenario of dark energy is investigated. The most popular model among them is 
■^J- . dubbed quintessence [9 . Besides, other scalar- field dark energy models have been studied, including phantom [6J, 
tachyon Q, quintom [8], ghost condensates @, etc. Also, there are other candidates, for example, Chaplygin gas 
which attempt to unify dark energy and dark matter [Toj . braneworld model (lTj and 5-dimensional gravity model 12 1 



which explain the acceleration through the assumption that spacetime has five dimensions instead of the usual four. In 
addition, since the cosmological scaling solution (i.e., the energy densities of dark energy and cold dark matter remain 
proportional) could probably alleviate the coincidence problem, interacting dark energy models are also proposed [l3[ . 

As we all know, observations at the level of the solar system severely constrain non-gravitational interactions of 
baryons, namely, non-minimal coupling between dark energy and ordinary matter fluids is strongly restricted by the 
experimental tests in the solar systems [l4j], we therefore neglect this possibility. However, since the nature of dark 
sectors remains unknown, it is possible to have non-gravitational interactions between dark energy and dark matter. 
So we focus on dark energy interacting with dark matter alone. 

Actually, many dark energy models are considered in the framework of classical Einstein cosmology. However, an 
outstanding problem in classical Einstein cosmology is the big bang singularity which is expected to be solved by 
quantum gravity. As a background independent quantization of general relativity, loop quantum gravity (LQG) is 
one of the best candidate theories of quantum gravity [TEj . It has been applied in cosmology to analyze our universe, 
known as Loop Quantum Cosmology (LQC) [lj|. In LQC, non-perturbative effects lead to —p 2 /p c corrections to the 
standard Friedmann equation and thus allow us the possibility of resolving any past and future singularities [H, [l?} ■ 
The modification becomes important when energy density of the universe becomes to be the same order of a critical 
density p c . When the correction term —p 2 /p c dominates during the evolution of our universe, it will cause the 
quantum bounce and hence avoid the singularity. Recently, more and more researchers have taken their attention to 
LQC for the appealing features: avoidance of various singularities [lH, inflation in LQC (T^, large scale effect |20| 
and so on. Concretely, some dark energy models are investigated in the background of LQC, such as phantom |21j . 
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coupling phantom [22j], quintom and hessence (23[, interacting dark energy model [Hj], etc. 

In this paper, we study the dynamical evolution of two classes of interacting dark energy models in classical 
Einstein and Loop Quantum Cosmology. Here some questions naturally arise as follows. Can these models alleviate 
the coincidence problem in classical Einstein cosmology? Are there scaling solutions arising from the effect of loop 
quantum cosmology? Can the future singularities be resolved in LQC? By our analysis, it turns out that in the former 
model, there are two attractors in classical Einstein cosmology and LQC. One is an accelerated scaling solution and 
the other is a baryon dominant solution. However, in the latter model, there exists one attractor in classical Einstein 
cosmology, which is a dark energy dominated solution rather than a scaling solution, whereas in LQC all fixed points 
are unstable. Thus, there exists no scaling solution in the latter case, namely, this kind of interacting dark energy 
model can not be regarded as a candidate to alleviate the coincidence problem. Also, we find that dynamical results 
in LQC are different from those in classical Einstein cosmology for both two kinds of interacting models. Our universe 
finally enters an oscillating phase in LQC. Moreover, the oscillating frequencies are significantly different for varied 
parameters of models. These results arc different from the those obtained in classical Einstein cosmology. Thus, LQC 
allow us the possibility of resolving future singularities. Hence, the quantum gravity effect may be manifested in large 
scale in the interacting dark energy models. 

In Sec. II, we study dynamical properties for the general case in classical Einstein cosmology and LQC. Then the 
dynamical results of two types of interacting models are respectively studied in Sees. Ill and IV. In Sec. V, the 
numerical results are presented. Finally, the conclusions are summarized in Sec. VI. 



II. INTERACTING DARK ENERGY MODEL IN CLASSICAL EINSTEIN COSMOLOGY AND LQC 



For a spatially flat universe, the total energy conservation equation is 

p + 3H(p + p)=0, (1) 

where H is the Hubble parameter, p is the total energy density and p is the total pressure of the background fluid. 

In our scenario, the universe contains dark energy, cold dark matter and baryonic matter. Moreover, the two dark 
sectors interact through the interaction term Q and the baryonic matter only interacts gravitationally with the dark 
sectors. Then the energy conservation equation is written as 

p b + 3Hp b = 0, (2) 



p m + 3Hp m = Q, (3) 

p d + 3H(1 + w d )p d = -Q, (4) 

where the subscripts b, m and d respectively denote baryonic matter, cold dark matter and dark energy. In this 
paper, we consider the simplest case of dark energy with constant equation of state uid = Pd/pd [25j . although the 
equation of state for dark energy could also be dynamic. Thus, Q denotes the energy density exchange in the dark 
sectors and the sign of Q determines the direction of energy transfer. A positive Q corresponds to the transfer of 
energy from dark energy to dark matter, while a negative Q represents the other way round. Due to the unknown 
nature of dark sectors, there is as yet no basis in fundamental theory for a special coupling between two dark sectors. 
So the interaction term Q discussed currently have to be chosen in a phenomenological way [2(jj . Since there is no 
clear consensus on the form of the coupling, different versions, that arise from a variety of motivations, coexist in the 
literature. 



A. Classical Einstein Cosmology 



In classical Einstein cosmology, the Friedmann equation is given by 

H 2 = —p= —{p d + Pm + Pb), (5) 

where k 2 = 8irG. Then differentiating the above equation with respect to cosmic time t and using the total energy 
conservation equation, we can get the Raychaudhuri equation 

H = -^{P + P) 

K 2 

= — 2"((1 + Wd)Pd + Pm + Pb)- (6) 
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To analyze the evolution of the dynamical system, we introduce the following set of dimensionless variables: 

= K 2 p d _ n 2 p m _ n 2 p b _ k 2 Q 

X ~ 3H 2 ' V ~ 3H 2 ' Z ~ 3H 2 ' U ~ 3ff 3 ■ 1 j 

Accordingly, the Friedmann constraint is 



and Eqs. (5) and (6) can be written as 



x + y + z = l, (8) 



-^2 = \{i+mx). (9) 

Furthermore, using these variables, the EoS of the total cosmic fluid is given by 

P w d% 

w = - = = WdX. (10) 

p x + y + z 

Then, inserting the expression (7) into Eqs. (2)- (6), we can obtain the following autonomous system: 

x = —?>Wdx(l — x) — u, (11) 
y = Swdxy + u, (12) 
z = Swdxz, (13) 

where the prime denotes a derivative with respect to N = In a. We set the current scale factor by ao = 1, then the 
current value of N reads Nq = 0. Setting x = y = z = 0, we can find the general solution of the critical points 
(x*, y», z*) of the autonomous system (11)-(13) as the type of (x* = 0, ii* = 0) and [x„ + y* = 1, z* = 0). 

B. Loop Quantum Cosmology 

Due to the quantum effects in LQC, we consider effective Friedmann equation with correction of the form [l6[ 

H 2 = ^p{l-^l (14) 
3 p c 

where p = pd + p m + Pb, the critical density p c = 16 ^? 3 Ppi measures the loop quantum effects, p p i is the Planck density, 
7 is the dimensionless Barbero-Immirzi parameter [27| . An important feature for the modified dynamics is that a p 2 
term which is relevant in the high energy regime is included in the classical Friedmann equation. The correction term 
predicts a bounce when the matter energy density reaches the critical value p c which is close to the Planck density. 
By the numerical simulation f]~6t | , it turns out that the modified Friedmann equation is valid in the whole evolutional 
trajectory of the universe including the bounce. Additionally, along with the total energy conservation equation, we 
get 

H = _£(p + „)(l- 2 -£-) 
2 p c 

= — —{{l+w d )Pd + pm + Pb){l-2 )• (15) 

2 p c 

Using the dimensionless variables defined in (7), the Friedmann constraint is 

, . . 3H 2 x + y + z , , „ 

(x + y + z)(l-— y - =1, 16 

n A p c 



and Eqs. (14) and (15) can be written as 



4 = !(2-*-y-*)(i + ^ 

H z 2 x + y + z 



, n (2-x-y-z){l+ _ , ° , J . (17) 
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Point 




Eigenvalues 




A 


( c 1 ~c 2 -w d -x s c 1 -c 2 +w d -x 3 n \ 


Ai = 3a; s , 

A2 = A3 = — |(ci — C2 — Wd — x s ) 


ei-e 2 — Uirf-Ks 


\ 2m d ' 2w d > U i 


2 


B 


( c 1 -c 2 — w d + x s c 1 —c 2 +w d + x s q\ 


Ai = — 

A2 = A3 = -|(C1 — C 2 — W d + Xs) 


ci-C2-« , d+ :c s 


1 2iv d ' 2w d > U ) 


2 


C 


(0,0,1) 


Ai =0, 

A2,3 = |(ci - c 2 - w d ± a; s ) 






TABLE I: The properties of the critical points for the interacting model I in classical Einstein cosmology. Here, the parameter 
x s is defined in Eq. (25). 



Furthermore, the EoS of the total cosmic fluid reads 

P w d x 

w = - = . (18) 

p x+y+z 

Then, Eqs.(2)-(4) combined with Eqs. (14)-(15) can be rewritten as the following autonomous system according to 
the expression (7), 

x = -3(1 + w d )x -u + 3x(2 -x-y- z)(l + - ™ d ^ + - ), (19) 

V = -3y + u + 3y(2-x-y-z)(l+ ™ dX ), (20) 

x + y + z 

z = -3z[l-{2-x-y-z)(l + ™ dX )]. (21) 

x + y + z 

The type of critical points of the autonomous system (19)-(21) can be summarized as (x* = 0,u* = 0), (2;* + y* — 
1,2* = 0) and (a;* ^ 0, y* = —(1 + w d )x*,u = -3(1 + w d )x*). 



III. INTERACTING DARK ENERGY MODEL I 



Let us first consider the interaction term Q = 3H(cip m + C2p d ) [28l 29], where c\ and Ci are coupling constants. 
Note that this form, which was first proposed in [3(| , is more general than those proposed in (2f| IU . The latter can 
be obtained from the former by setting c\ = c 2 = c, c\ — or c 2 = 0. According to Ref. 
constants c\ and c 2 have the same sign to achieve a physically viable model. 



281 ] . we assume the coupling 



A. Cosmological Dynamics in Classical Einstein Cosmology 

In classical Einstein cosmology, the autonomous system (11)-(13) can be written as 

x = -3[(l + c 2 +w d )x + c iy }+3x(l + w d x), (22) 
y = -3[(1 - cx)y - c 2 x] + 3y(l + w d x), (23) 
2 = 3w d xz. (24) 



Furthermore, we can obtain the critical points (a;*, y*, z*) of the autonomous system as follows: 

-2-W d -X s C 1 —C 2 +W a 

2w d ' 2w d 



. Point A: (- Ci -%-y- x % - Cl ~ C2 2 t?" a! - , 0), 



• Point B: f_ ci-c2-w d +x. _ c 1 -c 2+Wd +x Q) 



• Point C: 



(0, 0, 1). 
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Point 


(a;*, y*,z*) 


Eigenvalues 




A 


/ c 1 -c 2 -w d -x 3 c 1 -c 2 +iv d -x s n\ 


Ai = 3a: 3 , 

A2 = — f (Cl — C2 — lUd — £s)j 
A3 = -|(2 - Cl + C 2 + Wd + £„) 


ci-c 2 -m rf -a; s 


I 2w d ' 2w d ' U J 


2 


B 


( c 1 -c 2 -iu d +x a ci-c 2 +w d +x s p,\ 


Ai = — 3t 
A 2 = -|(ci - c 2 - w d + x 3 ), 
A3 = -|(2 - Cl + C 2 + Wd - x s ) 


Cl-c 2 -™ £ i + a; s 


V 2w d ' 2w d > U J 


2 


C 


(0,0,1) 


Ai = -3, 
A2,3 = |(ci — C2 — w d ± a; s ) 






TABLE II: The properties of the critical points for the interacting model I in LQC. Also, the parameter x s is defined in Eq. 
(25). 



Here the parameter x s is defined by 



= \J (ci - c 2 - w d ) 2 + Aw d c\. 



(25) 



To study the stability of the critical points for the autonomous system, we substitute linear perturbations x — > 
x* + Sx, y — > y* + <5y and z — > z» + <5z about the critical points into the autonomous system Eqs. (22)-(24). To 
first-order in the perturbations, we get the following evolution equations of the linear perturbations: 



Sx — — 3(c 2 + Wd — 2wd.x^)5x — 3ci5y, 
5y = 3(c 2 + w d y„)Sx + 3(ci + w d x*)5y, 
Sz = 3wdZ*5x + 3wdX*Sz. 



(26) 
(27) 
(28) 



The three eigenvalues of the coefficient matrix of Eqs. (26)-(28) determine the stability of the critical points. We list 
the three eigenvalues for each point in Table I. We examine the sign of the eigenvalues of points A and B and find 
that point A is not stable if it exists, whereas point B is stable if < c\ < —Wd and < c 2 < c± — w d — 2^—Wd.ci, i.e., 
the critical point B is always the stable attractor solution if it exists. In addition, at point B, from the expression of 

the acceleration condition w < — 1/3 shows that if 



the total EoS w* = - ^-^-w d +x 



w d > -1/3, 
-2/3 < w d < -1/3, 
w d < -2/3, 
w d < -1/3, 



< ci < —Wd, < c 2 < (1 + 'dw d )ci - w d 

w d + 2/3 < ci < -1/ (9w d ), < c 2 < (1 + 3w d )ci - w d -r ±/ ^ 
< ci < -l/(9wA < c 2 < (1 + 3w d )ci - w d + 1/3 

< c 2 < ci - w d - 2y/-W d Cl, 



J < ci < -l/(9io d ), 
- l/{9w d ) < ci < -w d 



1/3 
1/3 



(29) 



point B is an accelerated scaling attractor where the energy densities of dark energy and cold dark matter remain 
proportional. Thus it may alleviate the coincidence problem. The baryon dominated point C is stable if ci — c 2 — w d < 
and 4w d c\ < 0. This condition corresponds to ci > and c 2 > ci — Wd under the prior condition Wd < 0. 



B. Cosmological Dynamics in LQC 



We consider the autonomous system (19)-(21) in LQC. By inserting the concrete form of Q into Eqs.(19)-(21), the 
autonomous system can be expressed as 



x = -3[(1 + c 2 + w d )x + ay] + 3x(2 - x - y - z)(l 
y = -3[(l-c 1 )y-c 2 x]+3y(2-x-y-z)(l + 
z = -3z[l-(2-x-y-z)(l+ WdX )]. 



w d x 



x + y + z 
w d x 



x + y + z 



), 



x + y + z' 

The corresponding critical points (x*,y*,z*) of the autonomous system (30)-(32) are obtained as follows: 
. Point A: (- Cl ~ C2 2 '7' a: % - Cl - C ^ w w d d - X ' , 0), 



(30) 
(31) 
(32) 
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• Point B: r_ ci-c a - Wd +x _ Cl - Ca + Wd+x Q) 

2w d 2w d 



• Point C: (0, 0, 1). 

In order to study the stability of the critical points for the autonomous system (30)-(32), we obtain the following 
evolution equations of the linear perturbations: 



Wd,xl WdX*(2 — — y» — z*)(x* + 2y* + 22*) 1 

*c* y* -f~ z* (•£* y# ^*)^ 



oa; = -6[—l + c 2 + + ix* + y» + z* + - — — - — — — — — - — , 2 \ox 



2w d x\ 2u>dX* 

-3 ci + a;* + 7 -r \oy — 3x* 1 + ^ fc, (33) 

[x*+y* + z*y ' [x* + y* + z*y 



or ,, , i 2w d y* (y* + z») 

<ty = 3[c 2 - (1 + w d )y* + ■ ; -2 \dx 

[x* + y* + z*y 



-3 1 + ci - 1 + w d )x* - 2y., - z* + ■ ■ —}Sy - 3y* 1 + ■ ■ -^}Sz, (34) 

(x* + y* + z*) 2 [x*+y*+z*y 



o r-i | 2w d (y*+Z,Q 



g n , tw d x* 2w d x*{x* +y*) u , 

-32* 1 + ■ ■ -TrJdy - 3 -1 + (1 + w d )x* + y* + 2z« - ■ ■ r^\oz. (35) 

(a;*+y* + ^*) (a;*+y* + ^*) 

The corresponding three eigenvalues of the coefficient matrix of Eqs.(33)-(35) are listed in Table II. It is easy to examine 
that point A is not a stable point, while point B is stable if Wd > —1,0 < Cj < —10^,0 < c 2 < c\ — Wd — 2^J—w d c\ 
or Wd < —1,0 < ci < — l/w d , (1 + Wd){c\ — 1) < c 2 < Ci — — 2-^/— w^ci. We find that for the case of w > — 1, the 
stable regions of parameters in LQC are the same as those in classical Einstein cosmology. Additionally, at point B, 
the accelerated condition w < —1/3 can be expressed as 

( < ci < -w d , < c 2 < (1 + w d )c x -w d + 1/3 if w d > -1/3; 

w d + 2/3 < Cl <-l/(9w d ), < c 2 < (1 + io d )ca - w d + 1/3 if -2/3 < life < -1/3; 

< ci < -l/(Qw d ), < c 2 < (1 + 3w d ) Cl - w d + 1/3 if -1< w d < -2/3; 



-1/ (9iUd) < ci < < c 2 < ci - w d - 2 v /3 u^ci if -1 < w rf < -1/3; 

-l/(9w d ) < ci < -l/(3w d ), (1 + w d )(ci - 1) < c 2 < ci - w d - 2/3 if w d < -1; 

. < ci < -l/(9tu rf ), (l + u; d )(ci - 1) < c 2 < (1 + 3w d )ci - w d + 1/3 if w d < -1 



(36) 



so that point B is an accelerated scaling attractor for any case above, which provides a possibility to alleviate the 
coincidence problem. It should be noted that point C is stable for the same condition as the case in the classical 
Einstein cosmology, and therefore, we will neglect discussing this point in the section V. 

IV. INTERACTING DARK ENERGY MODEL II 

In this section, we study the interaction dark e nerg y model with a constant transfer rate: Q — Tp rn [32| . which 
has been used in reheating [33| , dark matter decay |34j curvaton decay [35| and the decay of superheavy dark matter 
particles to a quintessence scalar field (36j . 

In order to study the dynamical evolution, we additionally define a new dimcnsionlcss variable as 

« = (37) 
where Hq denotes the current value of the Hubble parameter, and for convenience, we introduce the parameter 

-no 
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Point 




Eigenvalues 




A 


(0,y,,l-y.,0) 


Ai = A 2 = 0, 
A3 = §, A4 = — 3utd 





B 


(1 01 


Ai = A 2 = A 3 = 3w d , 
\d = -d + tuwl 




C 


I 1 l+«>d n 3 \ 
V w d ' w d ' u > P> 


Ai = A2 = —3, 
A 3 ,4 = §(-l-w d ±\/«^-l) 


-1 



TABLE III: The properties of the critical points for the interacting model II in classical Einstein cosmology. 



A. Cosmological Dynamics in Classical Einstein Cosmology 



In classical Einstein cosmology, using Eqs.(2)-(6), (37) and (38), we concretely express the autonomous system 
(11)-(13) as 

x = —3wdx(l — x) — /3vy, (39) 

y = -3w d xy + /3vy, (40) 

z = 3wdxz, (41) 

' 3 

v = -v(l + w d x), (42) 

which has three critical points as follows: 
• Point A: (0, y*, 1 - j/*, 0), 



• Point B: (1, 0, 0, 0), 



• Point C: i±3ai, 0, -I). 

Substituting linear perturbations x — > + fe, y — > y» + Jy, z — > z* + <5z and w — > t>* + <5u about the critical 
points into the autonomous system Eqs.(39)-(42), to first-order in the perturbations, we get the following evolution 
equations of the linear perturbations: 

Sx — — 3wd(l — 2x*)8x — fiv^by — fiy*8v, (43) 

Sy = 3u>dy*5x + (3wd,x* + 0v*)5y + f3y*5v, (44) 

5z = 3wdz st 8x + 3wdX*5z, (45) 
, 3 3 

Sv = -WdV*Sx + -(1 + WdX*)5v. (46) 

The four eigenvalues of the coefficient matrix of the above equations determine the stability of the critical points. In 
Table III, we list the eigenvalues for each point. Then, it is clear that the critical points A and C are not stable if they 
exist, while point B is stable when u>d < — 1. Since the total EoS at point B is w» = u>d, point B is an accelerated 
attractor. Note that it is a dark energy dominated solution, rather than a scaling solution. Thus, this kind of 
interacting model in classical Einstein cosmology can not be regarded as a candidate to alleviate the coincidence 
problem. 
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Point 




Eigenvalues 


w* 


A 


(0,y.,l-i/.,0) 


Ai = 0, A 2 = -3 
A3 = f , A4 = — 3wd 





B 


(1,0,0,0) 


Ai = A 2 = 3w d , 
A 3 = |(1 + w d ), A 4 = -3(1 + w d ), 




C 


(a;.,-(l + w d )a;«,0, |) 


Ai = 0, A 2 = 3[Wd+2> 
A 3 ,4 = |{-1 - w d ± v/(l + tu d )[-3 + to d (l - 2a;,)]} 


-1 



TABLE IV: The properties of the critical points for the interacting model II in LQC. 



B. Cosmological Dynamics in LQC 



In LQC, inserting Eqs.(37)and (38) into the evolution equations (19)-(21) and making use of Eqs.(14) and (15), we 
get the autonomous system 

x = -3(l + w d )x- pvy + 3x(2-x-y- z)(l+ mX ), (47) 

x + y + z 

V = -3y + 0vy + 3y(2-x-y-z)(l+ ™ dX ), (48) 

x + y + z 

z = -3z[l -(2-x-y- z){l + ^f^)], (49) 

v = * vi 2- X -y- z)i l + -^-). (50) 

The critical points of the autonomous system (47)- (50) is obtained as 
• Point A: (0, y*, 1 - y*, 0), 



• Point B: (1, 0, 0, 0), 



• Point C: (a;*, -(1 + Wd)x*, 0, |). 



In order to study the stability of the critical points for the autonomous system (47)- (50), we obtain the following 
evolution equations of the linear perturbations: 

, ' or t , , o , , WdX * w d x*(2 - x* - y* - z»)(z* +2y* + 2zA 
<5.T = -3[-l + + 2x„ + y* + z* H — — }dx 



*y = -3y.[l + ^- ( ^ + ^ + ^ )2 ]^ 

+{/*«. + 3[1 - (1 + - 2y. - z, + 2 -^^±^]}S y 3y.[l + ^ ]fe + ft,.*;, (52) 

x ' q n 2w d (y» + z*) 

dz = -3z» 1 + w d - - da; 

(a;* + y* + z*y 



_3z4l + - l WdX * ]6y 3[-l + (1 + + y* + 2z» - ^g^±j^ ]fe, (53) 

(x» +y» + z*) 2 (x» + y»+z,) 2 



x ' 3 n , 2w d (y» + z„) 

*" = '2^ [1 + ^- (x, + y, + z,) 2 ]fa 

3 2w d a;, 3 2w d x st 3 w d x 

--«* 1 + 7 r^]Sy - -vJl + - ^]6z + -(2-x-y- z)(l H )Sv. (54) 

2 L (ar. + y* + z, 2J 2 L (ar. + y* + z, 2J 2 V y A x + y + z' v ; 



9 



Solving the four eigenvalues of the coefficient matrix of the above equations, we list them in Table IV. It is not difficult 
to see that the critical point A is not stable. For point C, to guarantee the energy densities of dark sectors to be 
positive, we get Wd < — 1 , and therefore the real parts of the eigenvalues A3 and A4 are positive. This means that 
point C is not stable. The critical point B is also not stable, since the sign of A3 is always opposite to the sign of A4. 
However, in classical Einstein cosmology point B is stable when Wd < —I. 

V. NUMERICAL RESULTS 

In what follows, we numerically study the dynamical results of the interacting dark energy models to confirm the 
complicated stability condition for the critical points in both interacting models. 

A. The Interacting Model I 

In the former model, there are two attractors in both classical Einstein cosmology and LQC. One is a baryon 
dominated solution which is neglected, and the other is an accelerated scaling solution. Numerical results for the 
interacting model I are shown in Figs. 1-8. 

In Fig.l, we depict the parameter space [c\, C2) to be stable by choosing Wd = —0.6 and Wd = —1.4. When Wd > — 1, 
the stable region in LQC is the same as that in classical Einstein cosmology, i.e., the critical point B is stable in the 
region I in both kinds of cosmology. However, when Wd < —I, point B is stable in the region I+II in classical Einstein 
cosmology, whereas in LQC, point B is an attractor only in the region II. 

In Fig. 2, we plot the phase space trajectories of the universe with Wd, c\ and c-i in the stable region. We find that 
the position of the critical point B depends on the EoS Wd and the coupling constants c\ and C2, but is independent 
of the theory describing our universe. However, the trajectories in the phase space depend not only upon Wd, c\ and 
C2, but also upon the selected theory. 

Fig. 3 shows the evolution of the total EoS w with the chosen parameters Wd, c\ and C2, satisfying the acceleration 
conditions (29) and (36). Apparently, we can see that in the final state the total EoS w tends to a constant, which 
depends on Wd, c\ and c 2 , but is independent of the theory describing our universe. 

Fig. 4 exhibits the trajectories of scalar factor a versus time for different values of parameters in LQC. We set k 2 = 1, 
and thus take p c = 1 since the value of p c is on the order of the Plank density, k -4 . From the figure, one can see that 
the evolution trajectories are significantly different for varied parameters. The bounce in scale factor occurs later for 
greater value of Wd, c\ or C2. Our universe finally enters an oscillating phase in LQC. 

In Figs. 5-8, we plot the evolution trajectories of the Hubble parameter H and energy density p versus time. The 
parameters we selected is in the unstable region in LQC. Differentiating Eq.(14) with respect to p, we find that H 
has a extremum value (dH/dp = 0) when p = p c /2. Additionally, the second order derivative of H reads 

(H > 0), (55) 
(H < 0). (56) 



Thus, H max — J K 2 p c /12 when H > 0, while H m i n = —^Jn 2 p c /12 when H < 0. In Figs. 5-8, with p c — 1, a calculation 
gives H max ^ min w ±0.2887 at p = p c /2. When H rts 0, we have the density p w p c and thus the bounce occurs. From 
the figures, we find that the expansion of our universe halts at the time when H w 0, and then contracts until H ss 
again. The universe goes on bouncing forward and backward. It is worthwhile to note that the oscillating frequencies 
of H{t) and p(t) depend upon the chosen coupling constants c\, c-i and EoS Wd- 

B. The Interacting Model II 

In the latter model, the dark energy dominated solution is the only attractor solution in classical Einstein cosmology, 
whereas there exists no attractor in LQC. Numerical results for the interacting model II are presented in Figs. 9-14. 
Since observations constrain the interaction to be sub-dominant today, which indicates |T| -C Hq, we select the 
parameter j3 to be very small in the following numerical analysis. 

In Fig. 9, we plot three-dimensional phase space trajectories of the universe with Wd = —1.2 in the stable region in 
classical Einstein cosmology. From the figure, we see that the trajectory curves from different initial conditions are 
converged at a point, the position of which is independent of any parameters. 
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In Fig. 10, the evolution of the total EoS w is plotted. It is easy to see that in the final state w tends to a constant, 
which equals to W4. It is worth noting that the evolution trajectory of w is not only independent of the theory 
describing our universe but also independent of the coupling constant (3. In other words, the coupling constant do 
not affect the evolution result. 

Fig. 11 show the trajectories of scalar factor a versus time for different values of parameters in LQC. We also take 
k 2 = 1 and p c = 1. The evolution trajectories for different values of parameters are distinct. The bounce in scale 
factor occurs later for greater value of Wd or (3. Our universe finally enters an oscillating phase in LQC. 

In Figs. 12-14, we plot the evolution trajectories of the Hubble parameter H and energy density p versus time. The 
parameters we selected as those above is in the unstable region in LQC. The universe goes on bouncing forward and 
backward. The oscillating frequencies of H{t) and p(t) depend not only upon EoS Wd but also upon the coupling 
constant /?. 

VI. CONCLUSIONS 

In previous sections, we have studied the cosmological evolution of two interacting dark energy models in classical 
Einstein and Loop Quantum Cosmology. We consider two kinds of interaction term between dark energy and cold 
dark matter. Note that observations at the level of the solar system severely constrain non-gravitational interactions 
of baryons. So the baryonic matter solely satisfies the energy conservation equation. By our analysis, we find that 
dynamical results in LQC are different from those in classical Einstein cosmology for both two kinds of interacting 
models. 

In the interacting model I, namely, Q = '&H(c\p m + c^pd), there are two attractors in both classical Einstein 
cosmology and LQC. One is a baryon dominated solution and the other is an accelerated scaling solution. Since the 
same results are obtained in both classical Einstein cosmology and LQC for the baryon dominated attractor, we only 
focus on the accelerated scaling solution. Interestingly, we find that if Wd > —1, the stable region in LQC is the same 
as that in classical Einstein cosmology, while when Wd < —1, the stable region in LQC is smaller than that in classical 
Einstein cosmology. The total EoS w approaches finally to a constant, which depends on EoS Wd and the coupling 
constants c\ and C2, but is independent of the theory describing our universe. When we select the parameters in the 
unstable region in LQC, the universe experience bouncing, which can resolve the singularity problem. The bounce 
in scale factor occurs later for greater value of Wd, c\ or C2. Furthermore, the oscillating frequencies are distinct for 
different parameters. 

However, in the interacting model II with Q = Tp m , there exists one attractor in classical Einstein cosmology for 
Wd < —1, which is a dark energy dominated solution rather than a scaling solution, whereas in LQC, all fixed points 
is not stable. Thus, there exists no scaling solutions in the interacting model II. So this kind of interacting dark 
energy model can not be regarded as a candidate to alleviate the coincidence problem in both classical Einstein and 
loop quantum cosmology. In classical Einstein cosmology, the final state is a constant, which equals to Wd and is 
independent of the coupling constant j3. The bounce in scale factor occurs later for greater value of Wd or j3. Our 
universe finally enters an oscillating phase in LQC. Moreover, the oscillating frequencies are significantly different for 
varied parameters. 

In summary, the interacting model I may alleviate the coincidence problem in both classical Einstein and loop 
quantum cosmology, depending on the values of the parameters selected in the model. However, the interacting 
model II can not be regarded as a candidate to alleviate the coincidence problem in both kinds of cosmology. Thus, 
dynamical results are different not only in different theories describing the universe but also in different interacting 
models. In addition, the results that our universe finally enters an oscillating phase in LQC, which are different 
from the those obtained in classical Einstein cosmology, show that LQC allow us the possibility of resolving future 
singularities. Therefore, the quantum gravity effect may be manifested in large scale in the interacting dark energy 
models. 
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FIG. 1: The stable regions in the (ci,C2) parameter space with the fixed Wd for the interacting model I. In the left plot 
(Wd > — 1), the critical point B is stable in the region I in classical Einstein cosmology and LQC. In the right plot (uid < — 1), 
point B is stable in the regions I+II in classical Einstein cosmology, whereas in LQC, point B is an attractor only in the region 
II. The region II in the left and the region III in the right represent the regions of the physically meaningless solution. 





FIG. 2: Three-dimensional phase space of (a;, y, z) with the fixed Wd for the interacting model I. The left and right two plots 
respectively denote classical Eintein cosmology and LQC. In the top two plots, we select the parameters ci =0.1 and C2 = 0.15. 
The bottom two plots is for ci = 0.1 and C2 = 0.5. 




FIG. 3: The evolution of the EoS of total cosmic fluid w with the fixed Wd for the interacting model I. The left is for classical 
Einstein cosmology and the right is for LQC. The parameter Ci is chosen as 0.1. 
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FIG. 4: The plots of scalar factor a as a function of time for the fixed parameters. The left plot is for Wd = —1.4 and ci = 0.1. 
The solid and dotted lines respectively correspond to C2 = 0.2 and C2 = 0.1. The middle plot is for Wd = —1.6 and C2 = 0.1. 
The solid and dotted lines respectively correspond to ci = 0.3 and ci =0.1. The right plot is for ci = C2 = 0.1. The solid and 
dotted lines respectively correspond to Wd = —1.4 and Wd = —1.6. 




FIG. 5: The evolution of the Hubble parameter H and the energy density p with respect to time with Wd = — 1.4, ci = 0.1, C2 = 
0.2. The solid, dashed, dash-dotted and dotted lines correspond to p, p d , p m and p b , respectively. 




FIG. 6: The evolution of the Hubble parameter H and the energy density p with respect to time with Wd = — 1.4, ci = 0.1, C2 = 
0.1. The solid, dashed, dash-dotted and dotted lines correspond to p, pd, p m and pb, respectively. 
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FIG. 7: The evolution of the Hubble parameter H and the energy density p with respect to time with Wd = — 1.6, ci = 0.1, C2 
0.1. The solid, dashed, dash-dotted and dotted lines correspond to p, pd, pm and pb, respectively. 
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FIG. 8: The evolution of the Hubble parameter H and the energy density p with respect to time with Wd = —1.6, c\ — 0.3, C2 
0.1. The solid, dashed, dash-dotted and dotted lines correspond to p, pd, pm and pb, respectively. 





FIG. 9: Three-dimensional phase space for the interacting model II with Wd = —1.2 in classical Einstein cosmology. The left 
is for the phase space of (x, y, z) and the right is for the phase space of (y, z, v). 
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FIG. 10: The evolution of the EoS of total cosmic fluid w for the interacting model II with j3 — 10 6 in LQC (the same as in 
classical Einstein cosmology). 



15 




FIG. 11: The plots of scalar factor a as a function of time for the interacting model II. The left plot is for Wd = —1.2, in which 
the solid and dotted lines denote the cases of /3 = 1CT 2 and /3 = 10 -6 , respectively. The right is for /3 = 10 -6 , in which the 
solid and dotted lines represent Wd = —1.2 and Wd = —1-4, respectively. 





16 



0.3 
0.2 
0.1 



-0.2 



w 



40 60 
t 



100 




FIG. 14: The evolution of the Hubble parameter H and the energy density p with respect to time for Wd 
The solid, dashed, dash-dotted and dotted lines correspond to p, pd, p m and pb, respectively. 
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